Low-energy monopole strength in exotic Nickel isotopes 
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Low-energy strength is predicted for the isoscalar monopole response of neutron-rich Ni isotopes, 
in calculations performed using the microscopic Skyrme HF+RPA and relativistic RHB+RQRPA 
models. Both models, although based on different energy density functionals, predict the occurrence 
of pronounced monopole states in the energy region between 10 MeV and 15 MeV, well separated 
from the isoscalar GMR. The analysis of transition densities and corresponding particle-hole con- 
figurations shows that these states represent almost pure neutron single hole-particle excitations. 
Even though their location is not modified with respect to the corresponding unperturbed states, 
their (Q)RPA strength is considerably enhanced by the residual interaction. The theoretical analysis 
predicts the gradual enhancement of low-energy monopole strength with neutron excess. 



The multipole response of nuclei far from the /3- 
stability line and the possible occurrence of exotic modes 
of excitation presents a growing field of research. Be- 
sides being intrinsically interesting as new structure phe- 
nomena, exotic modes of excitation might play an im- 
portant role in r- process nucleosynthesis [ij. Low- lying 
El strength has been measured in neutron-rich Oxygen, 
Neon, and Tin isotopes and more recently in 68 Ni 

p|. The interpretation of the dynamics of the observed 
low-energy El strength in nuclei with a pronounced neu- 
tron excess is currently under discussion (see Ref. [f| 
for a recent review). In light nuclei such as the Oxygen 
isotopes the onset of dipole strength in the low-energy re- 
gion is associated with non-resonant independent single- 
particle excitations of loosely bound neutrons. However, 
in the case of 68 Ni the low-lying dipole strength is found 
to be rather collective 0, and several theoretical anal- 
yses have predicted the existence of the pygmy dipole 
resonance (PDR) in medium-mass and heavy nuclei, i.e., 
the resonant oscillation of the weakly-bound neutron skin 
against the isospin saturated proton-neutron core. More 
recently it has been shown that the low-lying El strength 
can exhibit a non trivial pattern, separated into two seg- 
ments with different isospin character. The more pro- 
nounced pygmy structure at lower energy is composed of 
predominantly isoscalar states with surface-peaked tran- 
sition densities. At somewhat higher energy the calcu- 
lated El strength is primarily of isovector character, as 
expected for the low-energy tail of the giant dipole reso- 
nance [|[. 

Of course, not only the dipole, but also other mul- 
tipoles could exhibit pronounced low-energy strength in 
neutron-rich nuclei. A theoretical study of the occurrence 
of pygmy quadrupole resonances has recently been re- 
ported in the framework the quasiparticle-phonon model 
[9(, and the monopole response in very exotic nuclei 
such as 60 Ca was investigated using the self-consistent 
Hartree-Fock calculation plus the random phase approx- 
imation (RPA) with Skyrme interactions [l(|. It was 
shown that near the drip line the monopole response 



could develop a very pronounced structure at low energy. 
The present study considers more realistic examples of 
nuclei in which the low-energy monopole strength could 
be measured in the near future. For the Ni isotopes, in 
particular, calculations of the monopole response in 78 Ni 
with the Gogny-RPA did not predict pronounced low- 
lying strength [11], as evidenced from the percentage of 
the EWSR exhausted in the low-energy region. In Ref. 
fj~2j the monopole strength in Nickel isotopes was stud- 
ied using the Skyrme-QRPA approach with the SkM* 
functional. It is, however, difficult to discern a clear pat- 
tern in the low-energy region because of the large artifi- 
cial smoothing factor. Here we study in more detail the 
occurrence of low-lying isoscalar monopole strength in 
neutron-rich Ni isotopes. In 68 Ni low-lying El strength 
has recently been measured A new technique has 
been developed that enables measurement of monopole 
strength in unstable nuclei [l3j |. An experiment has 
very recently been performed at GANIL to measure the 
monopole strength in 68 Ni [14], and the analysis is in 
progress. 

To minimize the model dependence of results, both 
Skyrme and relativistic energy density functionals are 
employed in the present analysis. Skyrme Hartree-Fock 
+ RPA calculations are performed for the nuclei 68 Ni and 
78 Ni, for which pairing effects are expected to play a neg- 
ligible role in the monopole response. The description of 
Skyrme HF+RPA approach can be found in numerous 
references (l^| - [l7| and will not be detailed here. We only 
mention that in this work RPA equations are solved in 
configuration space, and the particle-hole space is chosen 
so that the energy-weighted sum rule is fully exhausted. 
The continuous part of the single-particle spectrum is 
discretized by diagonalizing the HF Hamiltonian in a 
harmonic oscillator basis. An important quantity that 
characterizes a given state v = (E V ,LJ) is its transition 
density: 
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FIG. 1. Skyrme-RPA isoscalar monopole strength functions 
in 68 Ni (in non-energy weighted sum rule units), calculated 
with the SLy4 (a) and SGII (b) energy density functionals 



with a corresponding definition of the neutron (proton) 
transition density 8p v n (8pp) when the summation in 
Eq. (p} is restricted to neutrons (protons). A smoothing 
factor of 600 keV is used in plots of strength distribu- 
tions. 

The monopole response of Ni isotopes is also ana- 
lyzed using the fully self-consistent relativistic quasiparti- 
cle random phase approximation (RQRPA) based on the 
Relativistic Hartree-Bogoliubov model (RHB) De- 
tails of the formalism can be found in Refs. [6l.ll8j. In the 
RHB+RQRPA model the effective interactions are im- 
plemented in a fully consistent way. In the particle-hole 
channel effective Lagrangians with density-dependent 
meson- nucleon couplings are employed [19j . and pairing 
correlations are described by the pairing part of the finite- 
range Gogny interaction. Both in the ph and pp channels, 
the same interactions are used in the RHB equations that 
determine the canonical quasiparticle basis, and in the 
matrix equations of the RQRPA. 

In Fig. [T]we display the isoscalar monopole strength 
in 68 Ni, calculated with the SLy4 (left) and SGII (right) 
non-relativistic functionals. The giant monopole reso- 
nance (GMR) is calculated at about 20 MeV, and both 
functionals predict a pronounced low-lying structure lo- 
cated between 13 and 16 MeV excitation energy. For 
the response calculated with SLy4, the proton and neu- 
tron transition densities of the three states that com- 
pose the low-energy monopole structure are plotted in 
Fig. [2j These states are not purely isoscalar: the tran- 
sition densities exhibit neutron dominated modes, and 
the proton and neutron densities are not in phase in 
the interior of the nucleus. The configuration analysis of 
these states shows that they correspond to almost pure 
single hole-particle excitations. A single configuration 
contributes with more than 98% to the total strength: 
neutron (2p 3 /2,3p3/ 2 ), (2pi/ 2 ,3pi/ 2 ), (lh/2,^5/2) for the 
states located at about 13.6 MeV, 14.6 MeV and 15.6 
MeV, respectively. In nuclei that do not exhibit a pro- 
nounced neutron excess these 2hu> unperturbed config- 
urations are located at higher energies, whereas in this 
case they are found below the GMR. The reason is that, 
because of the neutron excess in 68 Ni, the 3p2f states 
are located close to the Fermi level and, moreover, the 
2plf shell is calculated just below the Fermi level. This 
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FIG. 2. (Color online) Skyrme-RPA neutron (solid line) and 
proton (dashed line) transition densities for the three low- 
lying isoscalar monopole peaks at 13.6 MeV, 14.6 MeV and 
15.6 MeV in 68 Ni, calculated with the SLy4 functional. 



is not the case for the 3s2d shell. The 2Tluj sd configura- 
tions that build the GMR are not lowered in energy and, 
therefore, the giant resonance is located at higher energy. 

Fig. [3] shows the RQPRA prediction for the isoscalar 
monopole response in 68 Ni, calculated with the relativis- 
tic functional DD-ME2 [20] . Despite a small difference 
in the predicted position of the GMR, the results for the 
low-lying part are very similar to those obtained with 
the Skyrme functionals, and the configuration analysis 
leads to the same conclusion about the non-collective 
character of these low-lying excitations. The compari- 
son of the RQRPA spectrum with the unperturbed states 
shows that the residual interaction affects the strength, 
but not the location of the low-lying states. They are 
well separated from the GMR which, in this relatively 
light nucleus, is found at high energy. The correspond- 
ing proton and neutron transition densities are plotted in 
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FIG. 3. (Color online) Monopole response of Ni calculated 
using the RHB+RQRPA model with the DD-ME2 [H func- 
tional. In addition to the RQRPA discrete spectrum and the 
corresponding Lorentzian averaged curve, the unperturbed 
spectrum is also displayed. 



Fig. |U The transition densities for the low-lying states 
at 12.16 MeV, 13.38 MeV, and 15.42 MeV (these states 
have the same particle-hole structure as the correspond- 
ing ones in the Skyrme-RPA calculation) correspond to 
almost pure neutron modes, and the radial dependence is 
very different form the typical isoscalar GMR transition 
densities exhibited by the collective state at 18.94 MeV. 
The enhancement of the low-lying monopole strength in 
the RQRPA case compared to the RHB one is due to a 
small increase of the magnitude of the neutron transition 
density, the strength being a momentum of the transition 
density. Also, a slight admixture of neutron and proton 
configurations contribute to these states. 

Similar results are also obtained for 78 Ni (Fig. [SJ, 
for which both Skyrme-RPA and RHB+RQRPA calcu- 
lations predict pronounced strength in the energy region 
around 15 MeV excitation energy, well separated from 
the GMR. The low-lying strength is more pronounced 
than in the case of 68 Ni. The overall structure pre- 
dicted by the DD-ME2 functional is shifted to some- 
what lower energy compared to the one calculated with 
SLy4. The comparison between the unperturbed spec- 
trum and the full RQRPA response nicely illustrates how 
the residual interaction builds the collective GMR from 
the states above 15 MeV excitation energy. The two low- 
lying unperturbed states are not lowered in energy, even 
though their strength is considerably increased by the in- 
clusion of the residual interaction in the full relativistic 
RPA. To illustrate the evolution of low-energy monopole 
strength in Ni isotopes, in Fig. [6] we show the monopole 
response of even-A 60_78 Ni isotopes calculated with the 
RHB+RQRPA. One can clearly follow how the low-lying 
strength in the energy region between 10 MeV and 15 
MeV develops with the neutron excess. Since these 
states, predicted both by Skyrme-RPA and RQRPA cal- 
culations, are non-collective, their occurrence essentially 
depends on the position of the neutron Fermi level in the 
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FIG. 4. (Color online) 68 Ni proton and neutron transition 
densities of the three low-lying RHB+RQRPA states at 12.16 
MeV, 13.38 MeV, and 15.42 MeV, compared to those of the 
GMR at 18.94 MeV. 



single-neutron spectrum. The RPA strength contained 
in these states, however, is markedly enhanced by the 
residual interaction. These predictions point to a possi- 
bly very interesting measurement of low-lying isoscalar 
monopole states in neutron-rich Ni isotopes, that would 
probe almost pure single-neutron configurations. A reso- 
lution of approximately 2 MeV can currently be achieved 
by the experimental setup dedicated to the measurement 
of the monopole response in unstable nuclei [13j . An 
improvement to a 1 MeV energy resolution is expected 
[lj], and the developpement of next generation active 
targets such as ACTAR [5l[ will proceed along this path. 
It should be noted, however, that the physical width of 
the monopole response (the spreading width, the decay 
width, and the Landau damping) could also prevent a 
clear separation between different low-energy states. 

In conclusion, fully self-consistent microscopic Skyrme 
HF+RPA and relativistic RHB+RQRPA calculations 
have been performed for the isoscalar monopole response 
in neutron-rich Ni isotopes, some of which should become 
experimentally accessible in the near future. Both the 
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FIG. 5. (Color online) Monopole response of 78 Ni calcu- 
lated using the Skyrme-RPA model with the SLy4 func- 
tional (in non-energy weighted sum rule units) (a), and the 
RHB+RQRPA model with the DD-ME2 functional (b). 
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non-relativistic and relativistic models, based on Skyrme 
functionals and the relativistic functional DD-ME2, re- 
spectively, predict the occurrence of pronounced low- 
energy monopole states in the energy region between 
10 MeV and 15 MeV, well separated from the isoscalar 
GMR. From the analysis of transition densities, transi- 
tion matrix elements and the particle-hole configurations 
that correspond to these states, it is evident that the 
low-energy monopole states represent almost pure neu- 
tron single hole-particle excitations. Their occurrence 
can be related to the closeness of the corresponding or- 
bitals to the neutron Fermi level in nuclei with large neu- 
tron excess. Even though their location is not modified 
with respect to the corresponding unperturbed states, 
their (Q)RPA strength is considerably enhanced by the 
residual interaction. The theoretical analysis predicts an 
increase of low-energy monopole strength with neutron 
excess. A measurement of these states would provide a 
direct information on the 2huj gap, and directly probe 
the single-particle spectrum in exotic neutron-rich nu- 
clei, that is known to be driven by the spin-orbit and 
the tensor terms of the effective inter-nucleon interaction. 
In particular, measurements of the monopole strength 
in neutron-rich Nickel isotopes are currently being per- 
formed, and this work will help to interpret the results. 
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